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Abstract: We study the hydro dynamic properties of strongly coupled SU{N) 
Yang-Mills theory of the Dl-brane at finite temperature and at a non-zero density 
of R-charge in the framework of gauge/gravity duality. The gravity dual description 
involves a charged black hole solution of an Einstein-Maxwell-dilaton system in 3 
dimensions which is obtained by a consistent truncation of the spinning Dl-brane in 
10 dimensions. We evaluate thermal and electrical conductivity as well as the bulk 
viscosity as a function of the chemical potential conjugate to the R-charges of the 
Dl-brane. We show that the ratio of bulk viscosity to entropy density is indepen- 
dent of the chemical potential and is equal to I/Att. The thermal conductivity and 
bulk viscosity obey a relationship similar to the Wiedemann-Franz law. We show 
that at the boundary of thermodynamic stability, the charge diffusion mode becomes 
unstable and the transport coefficients exhibit critical behaviour. Our method for 
evaluating the transport coefficients relies on expressing the second order differential 
equations in terms of a first order equation which dictates the radial evolution of the 
transport coefficient. The radial evolution equations can be solved exactly for the 
transport coefficients of our interest. We observe that transport coefficients of the 
Dl-brane theory are related to that of the M2-brane by an overall proportionality 
constant which sets the dimensions. 
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1. Introduction 



There has been recent interest in constructing holographic duals which model phe- 
nomena and properties observed in macroscopic low energy physics. Such holographic 
duals may provide new insights because the properties and phenomena of interest 
usually lie in a regime which is strongly coupled in the field theory description but 
semi-classical from the gravity point of view. Transport properties of various systems 
which admit holographic duals have been evaluated from the gravity description. A 
universal result which has emerged out of these investigations is that the ratio of 
shear viscosity r] to the entropy density s for field theories which admit gravity duals 
in the two derivative approximation is given by [1,2] ^ 

" ' (1.1) 



S 4:7lk 



B 

where h is the Planck's constant and ks is the Boltzmann's constant. This ratio has 
been evaluated for well known AdS/CFT pairs like A/" = 4 super Yang-Mills as well as 
simple phenomenological gravity models. Other gauge/gravity duals which involve 
near horizon geometries which are not asymptotically anti-de Sitter backgrounds 
like that of Dp-branes, p > 2 [4, 5] have also been studied. This ratio for these 
backgrounds has also been shown to be h/AirkB [6]. 

In [7], we began an investigation of macroscopic properties of the 1 + 1 dimen- 
sional field theory of the Dl-branes. In 1 + 1 dimensions, there is no shear, therefore 
it is necessary to study non-conformal field theories to obtain non-trivial hydrody- 
namic coefficients. Dl-branes are interesting as they provide the simplest and the 
most symmetric non-conformal 1 + 1 dimensional field theory which admits a gravity 
dual. The theory is the maximally supersymmetric Yang-Mills with SU{N) gauge 
group. It can be obtained as a dimensional reduction of A/" = 4 SYM from 3 + 1 
dimensions. In [7], we isolated the sound mode in gravity and evaluated the speed 
of sound Vs and the bulk viscosity ( in the following regimes 

(i) < T < v^, (1.2) 

(a) VxN^^ < T < yXiv-^/l 

Here, A = Qym^ is the t' Hooft coupling and T is the temperature. In the above 
regimes, the field theory of the Dl-branes admits a gravity dual [5] which for the 
purposes of evaluating transport coefficients reduces to an Einstein-dilaton theory in 
3 dimensions. In [7], it was shown that 



^See [3] for a recent review and list of references on related topics. 
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for hydrodynamics of the Dl-brane theory, here c is the velocity of hght. It was also 
seen that theories arising form Dl-branes at cones over Sasaki- Einstein 7- manifolds 
give rise to the values in (1.3). It was suggested that there might be a class of non- 
conformal field theories which admit 3d gravitational backgrounds for which C/s = 
h/AnkB- For the rest of the paper, we will work with units in which h = = c = 1. 

In [8], it was shown that the supergravity fluctuations which determine the hy- 
drodynamic coefficients of the uncharged Dl-branes were related by dimensional 
reduction to that of the M2-branes. The dimensional reduction related the shear 
viscosity of the conformal hydrodynamics of the M2-brane to that of the bulk vis- 
cosity of Dl-branes. This explained why C,/s = l/47r, it can essentially be traced to 
the relation (1.1) for the M2-branes. It also explained why the value of the speed of 
sound of the Dl-brane theory behaves as though it is a conformal theory in 2 -|- 1 di- 
mensions. One expects a similar connection for the transport coefficients between the 
Dl-brane theory with finite charge density and the corresponding M2-brane theory. 
This would imply that the ratio C,/s will be independent of the chemical potential 
and continues to be l/47r since it is related to the ratio rj/s oi the M2-brane theory. 
There should also be similar relationships between other transport coefficients like 
conductivity. This is one of our motivations to explore the hydrodynamics of charged 
Dl-branes. There is a need to develop novel theories for 1-1-1 dimensional condensed 
matter systems as many higher dimensional models can't be applied here and there 
is a profusion of knowledge through experiments about new such systems and their 
properties [9]. So another reason is to study the macroscopic properties of strongly 
coupled 1-1-1 dimensional field theories which admit gravitational duals. Gravity 
duals of 1 -|- 1 dimensional systems with a well defined field theory have not been 
extensively studied ^. These systems play an important role in many quantum phe- 
nomena and it is worthwhile to see what insights the gauge/gravity correspondence 
gives in this context with a well defined field theory in mind. 

In this paper, we study the hydrodynamics of Dl-branes at finite charge density 
in a regime which admits a gravity description. The gravity dual description involves 
a charged black hole in an Einstein-dilaton-Maxwell scalar system in 3 dimensions 
which is obtained by a consistent truncation of spinning Dl-branes in 10 dimensions. 
We study two situations: 

1. The case in which the charge density corresponding to a single U{1) of the 
5*0(8) R-symmetry of the Dl-brane theory is turned on, we call this the single 
charged Dl-brane. 

2. The situation in which equal charge densities along the 4 Cartan's are turned 
on, we call this the equal charged Dl-brane. 

^Holographic duals of 1+1 dimensional systems from a bottom up approach without a known 
boundary field theory were studied in [10, 11]. 
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In both these cases, we see that both the speed of sound and the ratio of bulk viscosity 
to entropy density is given by (1.3). The values of these quantities are independent of 
the chemical potential. We also evaluate the charge conductivity, the charge diffusion 
constant, the sound diffusion constant and the thermal conductivity for both the 
situations and compare the results for which the corresponding M2-brane calculation 
has been done. We see that apart from an overall proportionality constant which 
sets the dimensionality of the transport coefficients in these theories, the transport 
coefficients are identical in the two theories. The results are summarized in the 
following table. 



Transport 
Coefficients 


Single-charged 
Dl brane 


Equal-charged 
Dl brane 


Equal-charged 
M2 brane 




1 


1 {.o—ky 


1 K^-ky 


167rG3 9v^T+fc 


167rG3 9(l+/c)2 


167rG4 9(l+fc)2 


C 


167rG3l,4 V 1 + k 


167rG3L4(l + ^) 




V 






r|,(l+fc)2 
167rG4-L'-* 


Dc 


L3(3-2fe) 


L3(fc+3) 
6r2j(l+fc)2 














12r|j(l+fc)2 


12r|j(l+fc)2 




rjj {2fc+3)(l+fc) 


r\ (3-fc)(l+fe) 


t\ (3-fc)(l+fc) 


8LG3 k 


8LG3 k 


8L'G4 k 



Table 1. Transport coefficients of Dl-branes and M2-branes. 

Til : radius of the horizon k: (R-charge)^ in units of Th- 

G3, G4: Newton's constant in 3 and 4 dimensions. 

L,L': radius of the orthogonal S"^ for Dl, M2-branes. 

(Jdc: electrical conductivity, (: bulk viscosity, 77: shear viscosity. 

DcDg-. charge diffusivity, sound diffusivity, kj-: thermal conductivity. 

Hydrodynamics of uncharged M2-branes were first studied in [12, 13]. We ob- 
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tained the shear viscosity of the charged M2-branes from the fact that i]/ s = I/Att 
[14]. The conductivity of charged M2-branes was obtained from [15]. The charge 
diffusion constant for the M2-branes at non-zero chemical potential has not yet been 
evaluated in the literature as far as we are aware. However for M2-branes at zero 
chemical potential, the charge diffusion constant has been evaluated in [16] ^ and 
it agrees with the k = limit of the Dl-brane theory answer. The sound diffusion 
constant for the charged M2-branes has been calculated by using Dg = ri/2{e + p) 
where e,p are the energy density and the pressure of the M2-branes. Notice that 
the bulk viscosity of the Dl-brane theory is proportional to the shear viscosity of 
the M2-brane theory. Another observation of our study of the transport coefficients 
of the charged Dl-brane is the following relationship between the bulk viscosity and 
the thermal conductivity 

^ = 4vr^ (1.4) 

where kt is the thermal conductivity, T the temperature and fi the chemical po- 
tential. This relationship is analogous to the Wiedemann-Franz law and a similar 
relationship between shear viscosity and the thermal conductivity has been observed 
in the case of single charged D3-branes [17]. Since the charged Dl-brane theory is ob- 
tained as a consistent truncation of spinning Dl-branes, there is a maximum allowed 
spin or charge beyond which the solution is thermodynamically unstable [18]. We 
show that the transport coefficients exhibit critical behaviour at the boundary of the 
thermodynamical instability. For the single charged case, we observe that the charge 
diffusion mode becomes unstable at the boundary of instability. This suggests that 
for this case, the thermodynamical instability can be better understood by studying 
the charged diffusion mode in more detail. 

This paper is organized as follows: In the next section, we introduce the single 
charged Dl-brane background and obtain the consistent truncation of the solution 
to 3 dimensions. We also review the thermodynamics of this solution and obtain 
the boundary of thermodynamic instability. In section 3, we study hydrodynamics 
of a charged fluid in 1 -|- 1 dimensions and obtain the dispersion relations of the 
two hydrodynamic modes, the charge diffusion mode and the sound mode in terms 
of thermodynamic variables. We then use the thermodynamics of the Dl-brane 
solution to explicitly evaluate the dispersion relations. We also determine the form 
of the retarded correlation functions of the stress tensor and the charge current 
using conservation laws. In section 4, we study the supergravity fluctuations of 
the single charged Dl-brane solution and isolate the gauge invariant fluctuations 
which correspond to the two hydrodynamical modes in the field theory. In section 
5, we determine the various transport coefficients from gravity using the relevant 
Kubo's formula. To do this, we reduce the problem to solving a set of coupled first 



^See below equation (3.32) in [16]. 
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order non-linear differential equations which are exactly solvable in limit required 
by the Kubo's formulae. These first order equations dictate the radial evolution of 
the transport coefficient. In section 6, we discuss the properties of the transport 
coefficients, their behaviour at the boundary of thermodynamic instability. We then 
discuss the connection of the Dl-brane theory to that of M2-branes. It will be 
interesting to compare our results with what is known for these systems. Appendix 
A contains the details of the consistent truncation which is required to obtain the 
charged Dl-brane solution in 3 dimensions. Appendix B contains the details of the 
evaluation of the transport coefficients for the equal charged Dl-brane. 



2. The R-charged Dl-brane 

In this section, we introduce the gravity dual of SU{N) Yang-Mills with 16 super- 
charges in 1 + 1 dimensions at finite R-charge density and state its domain of validity. 
We then discuss its thermodynamic properties. This section will also serve to set up 
notations and conventions. 

In [5], it was argued that SU{N) Yang-Mills with 16 supercharges in 1 -|- 1 
dimensions at large is dual to the near horizon geometry of A^ Dl-branes. Heating 
up this theory to a finite temperature T, the gravity dual is given in terms of the 
near horizon geometry of non-extremal Dl-brane solution. The gravity dual can be 
trusted in the domain 

v^A^~i « T « v^, (2.1) 



where A = \/ Qym^ t'Hooft coupling of the theory. The only non-trivial 

viscous transport coefficient of this system was evaluated using this gravitational 
dual in [7]. We now wish to turn on finite i?-charge density in the field theory. By 
the usual gauge/gravity correspondence, the 5*0(8) isometry of the S'^ present in the 
near horizon geometry of the Dl-branes corresponds to the 5'0(8) R-symmetry of 
the Yang-Mills. Therefore to turn on i?-charge density, it is necessary to consider 
Dl-branes with angular momentum. The near horizon supergravity solution of non- 
extremal Dl-branes spinning along one of the Cartan directions of 50(8) is given by 
[18]. 

ds^ = Hf'\-fde + dz"") - 2Hf^^lsm^ Odtdct), 

+hI'^ [Idr^ + r2(Arfe2 + H sin^ ed<p' + cos' Od^l 
= Hl'\ 

^(2) sin2 ed<j^ A dz, (2.2) 

Hi L'^ 



-6 - 



where 



/2 cos^ e f 
A = l + ^, H=l + l„ (2.3) 



Ar6 Ar6 



1 



The above solution is written in the Einstein frame. dQf is the metric of a unit 
5-sphere and 

= gU2VNia')\ = ^ (2-4) 

with gs, a' being the string couphng and the string length respectively. A^'^^ is the 
gauge potential for the RR 2-form sourced by the Dl-branes. Note that the above 
solution reduces to the non-spinning near horizon solution of the non-extremal Dl- 
brane when one sets the angular velocity I = 0. For completeness, we mention that 
the background in (2.2) is a solution of type IIB supergravity equations of motion in 
10 dimensions obtained from the following action 



1 ' .10. 



R-ldM<pd'U-^e^F,f 



(2.5) 



To study the hydrodynamics of this solution, one needs to consider perturba- 
tions of this solution along the brane directions {t, z) and the radial direction. The 
fluctuations along the 7-sphere do not play any role. Thus to simplify our analysis, 
it is convenient to perform a Kaluza-Klein reduction of this solution to 3 dimen- 
sions. Using the results of [19], it can be shown that the 10 dimensional solution in 
(2.2) admits a consistent reduction on the S"^ sphere to the following solution in 3 
dimensions 

ds^ = [-cldt^ + c^dz"^ + 4c/r2) , (2.6) 



A = -T^Z^^> = -31og(-), VP 



Here H and K are defined as 

72 /2 6 

The details of this Kaluza-Klein reduction are given in Appendix A. The rotation 
along one of the Cartan directions reduces to the charge denoted by the gauge po- 
tential At in 3 dimensions. Note that the deformation of round 5^ metric in (2.2) 
parametrized by A results in an additional scalar \1/ in 3 dimensions. It can also be 
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shown using this consistent reduction that the background in (2.6) is a solution of 
the equations of motion of the following action 

-^^^^^'^ + ^^M*^^"^ + ^^^'(1 + ^"')) ' (2-8) 

where 

k = t^' = ^"-"aliaf. (2.9) 

Thus the 10 dimensional rotating Dl-brane solution reduces to a charged black hole 
of an Einstein-Maxwell-dilaton system along with a scalar. The R-charge is given by 
the gauge potential At corresponds to rotation along the S*^ in 10 dimensions. As a 
simple consistency check, note that both the action in (2.8) and the solution in (2.6) 
reduces to the truncation studied in [7] ^ for the uncharged Dl-brane. Since the above 
solution is a consistent truncation to 3 dimensions, any solution to hydrodynamic 
fluctuations studied in 3 dimensions can by lifted to 10 dimensions. For completeness, 
we write down the equations of motion of the action given in (2.8). 



G ^Cj^jy^ + C 0, 

+ l^e^<^/3(i + vl/-i), 



1 1 



□0+4e'^/3(2 + *"') = O, 

□ logvl/ - "^e-^^I^F^^F^^ + ^^e^^l\\ - 



2.10) 
2.11) 
2.12) 
2.13) 
2.14) 

2.15) 
2.16) 



We refer to the solution in (2.6) as the single charged Dl-brane. The equal charged 
Dl-brane solution in which equal charge density along all the 4 Cartans of the 5*0(8) 
are turned on is given in (A. 20) of Appendix A. 

2.1 Thermodynamics of the R-charged branes 

The thermodynamic properties of spinning D-branes were studied in complete gen- 
erality in [18] from which we can read out the thermodynamic properties of the 



*See equations (4.3), (4.5), (4.6), (4.7). 
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black hole of interest given in (2.6). We now summarize the relevant thermodynamic 
properties. The Hawking temperature and the entropy density are given by 

where k is given by 



(2.18) 

and th is the radius of the horizon which is given by the largest root of the equation 

r% + r%l^-4 = 0. (2.19) 
The energy density and the free energy density is given by 



e 



Here we have also identified the pressure using its thermodynamic relationship with 
free energy density. The charge density p and its conjugate the chemical potential /i 
are given by 

P = 8^^' ^ = Mr)\r^oo - A{r)\r, = (2.21) 

Note that we have defined the chemical potential as the voltage difference between 
the boundary r — )■ oo and the horizon. In writing these thermodynamic quantities, 
we have used the relation (2.19). 

For the black hole solution given in (2.6) with very large charge, there exists a 
thermodynamic instability. This instability is equivalent to the instability occurring 
in Dl-branes which are rotating too fast [18]. Given the energy density of the system, 
the thermodynamic stability is determined by the condition 

Hs = det ^ > 0. (2.22) 

\ dpds dp'^ / 

To evaluate it, it is convenient to write the above Hessian as 

\drQ dl dl dro J \drQ dl dl dr^ J 

where we have used the chain rule and standard thermodynamic relations. Using 
the expressions for the thermodynamic variables given in (2.17), ( 2.20) and (2.21), 
it can be shown that the Hessian reduces to 
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Thus the condition for thermodynamic stabihty imphes the following restriction on 
the values of the R charge 

k < I (2.25) 

Finally, for completeness, we mention that the condition for the validity of the su- 
pergravity solution of the non-extremal spinning Dl-brane remains the same as that 
of the non-extremal brane and is given by 



^N--^ « T « V\. (2.26) 
The bound k < 3/2 in terms of field theory chemical potential can be written as 

Therefore the transport coefficients evaluated in this paper are valid in the regime 
given by (2.26) and (2.27) of the field theory. 



3. Hydrodynamics of a charged fluid in 1 + 1 dimensions 

In this section, we show that a charged fluid in 1 + 1 dimensions has two hydrodynamic 
modes and derive their dispersion relation. The stress tensor and the charge current 
of a relativistic fluid in 1 + 1 dimensions are given by 

T>"' = (e + p)u>'u'' + Pr]"" - C(m^m" + v'"ldxu\ (3.1) 

where u'^ is the 2- velocity with m^m'^ = — 1 and ( is the bulk viscosity and a the 
conductivity. The remaining variables e,p,p,fi refer to the energy density, pressure, 
charge density and the chemical potential of the system respectively, t]'^'^ refers to 
the flat Minkowski metric with signature (—1, 1). The equations of motion of the 
fluid are given by the following conservation laws 

d^T'^" = 0, d^r = 0. (3.2) 

We now wish to obtain the linearized hydrodynamics modes, therefore let us consider 
small fluctuations from the rest frame of the fluid. The 2-velocity is then given by 

n° = 1, = 5u'. (3.3) 

Note that m° = 1 up to the linear order due to the constraint m'^m^ = — 1. In 
considering the small fluctuations, one should keep in mind that spatial and temporal 
variations of the thermodynamic quantities are all of linear order. We can write the 
stress energy tensor to the linear order as given below 

T°° = e + (5T°°, T^z^^r^oz^ ST'' = p —dJT^'. (3.4) 

e+p 
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In writing the above form of the stress tensor, we have ehminated 5u^ using 

5u' = — , dM' = (3.5) 
e+p e+p 

As we are working only to the hnear order on taking the spatial derivative of 6ux, the 
derivative acts only on 6T^^. This is because derivatives of thermodynamic quantities 
are first order and therefore contribute only at second order in the above equation. 
Similarly the current density can be written as 

f = P + Sf, f = 6f = p— - aTdzfi, (3.6) 

e + p 

where ft = /i/T and we have again used (3.5). It is convenient to work with ther- 
modynamic variables in which the energy density e and the charge density p are the 
independent variables and all other thermodynamic quantities are functions of e and 
p. Then we can write Sj^ as 

Sf = p— - aT {d^pdzST^^^ + dppdzSf) . (3.7) 
e + p ^ ' 

Substituting the form of the stress tensor and the current density given in (3.4), (3.6) 
and ( 3.7) into the conservation equations (3.2), we obtain 

5o5j° + - oT {d^pdl 5T°° + d,-pdlb3') = 0, (3.8) 

ao5T°° + a,<5T°" = 0, 



doST""' + 7f ^^^T"" + ^djf ^diST^' = 0. 

\de op J e + p 

The above three equations determine the linearised hydrodynamic modes. Perform- 
ing the Fourier transform of the equations given in (3.8) both in position and time, 
we obtain the following set of algebraic equations 

{-zoo + aTdpM^)6f + ^<^^°' + crTd^M^ST^^ = 0, (3.9) 

-2u;5T°° + iq5T^' = 0, 
iqd,p6T^^ + iqdppSf + i-iuj + ) 5T°^ = 0. 

The above equations have non-trivial solutions for the fluctuations 5j^,5T^^,5T^^ 
only if the following constraint on u is satisfied. 



-lu + aTq%p) tu' - q'd.p + + q'dpP ^— + aTq'd.p = 0. 

' e + p J \e + p ' 



(3.10) 
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To solve for uj in terms of g, we can assume the following expansions for u 



0J = Vsq — iDsq +■■■, oo = —iDcq +■••. (S-H) 

Substituting the first expansion of u in terms of q given in the above equation in the 
constraint (3.10) and matching terms of 0{q^) and O(g^), we obtain the following 
expressions for the sound speed and its damping coefficient 

vl = (8.p+-^a,p), (3.12) 

Similarly substituting the second expansion for u given in (3.11) in the constraint 
(3.10) and demanding that the leading coefficient of O(g^) vanishes, we obtain the 
following value for the charge diffusion constant Dc 

= ar ^-^f-^-y^ . (3.14) 



deP + p 



HP 



It can be shown that these are the only two modes of the equations of motion of 
linearized hydrodynamics. To summarize, the two modes are the sound mode and 
the charge diffusion mode given by the dispersion relations in ( 3.11). 

We can now use the thermodynamic properties of the charged black hole given 
in (2.17), (2.20) and (2.21) to evaluate the dispersion relations explicitly. From 
(2.20), note that the pressure just depends on the free energy of the system and 
is independent of the charge density. Therefore for the R-charged Dl-brane, the 
dispersion relations simplify to 



00 = —zai — 
op 



We can further simplify the charge diffusion constant as follows 



D, = a[dppi-^dpT), (3.16) 



T 

dip. _ p. diT 

3L3 (3-2A;) 



cr(167rG 



3 



2rl {3 + 2ky 



To obtain the second line, we have used chain rule and also the fact that the en- 
ergy density e is independent of /. The last line is obtained by evaluating all the 
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the derivatives of the thermodynamic quantities using ( 2.17), ( 2.20) and ( 2.21). 
Therefore we see that the charge diffusion mode is given by 

(3.17, 

Note that if the conductivity a remains finite at the boundary of thermodynamic 
stabihty k = 3/2, the charge diffusion mode becomes unstable. Later in this paper 
we will explicitly evaluate the conductivity of the charged Dl-brane solution and 
show that it is indeed finite at = 3/2 and thus at the boundary of thermodynamic 
stability, the charge diffusion mode becomes unstable. 

One way of reading out the transport coefficients is to study the hydrodynamic 
modes and identify the coefficient of the dissipative parts. From (3.15) and (3.17), 
we see that we can read out both the bulk viscosity and the conductivity. Another 
approach is to use Kubo's formula which directly give the transport coefficients in 
terms of the two point functions. Let us first define the various retarded Green's 
functions: 

G^.ap{co.q) = -i j rf2x^(t)e-^('^*+''^)([T^,(x),r,^(0)]), (3.18) 
G^,p(a;,g) = -i j ci2x^(t)e-^('^*+^^)([J^(x), T,,(0)]), 

G,,{uj,q) = -i j <fxe{t)e-'^^'^'^'\[J,{x),JM])- 
Conservation laws and symmetries constrain the form of G ^papi'^ ■, q) to be [7] 

G^ua^iu, q) = P^^uPai^Gsiuj, q) , (3.19) 

where P^i^ is defined by 

P,u = - (3.20) 

and = {—u, q). Thus the two point function of the stress tensor is determined just 
by one function Gb- For future reference, we write down the following component of 
this correlator 

Gzzzz = :r^GB{uJ,q)- (3.21) 

Similarly one can show that conservation laws k^G^^{uj,q) = determine the form 
of the retarded two point function of the currents to be [20] 

G^,{co,q) = P^,Gj{co,q). (3.22) 

We write down the following component of this two point function 

GU^,q) = ^^Gjicu,q). (3.23) 
io^ — q^ 
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What is left now is the retarded two point function of the stress tensor and the charge 
current. Though we will not be requiring the form of this two point function, for 
completeness, we state that conservation laws and symmetries determine this two 
point function to be 

G^,p{uj, q) = e^^k'^P^.GsiuJ, q), (3.24) 

where e^j, is the antisymmetric tensor with etz = —tzt = 1- 

The transport coefficients, bulk viscosity C and the conductivity a are given by 
the following Kubo's formulae 

C = lim g = 0) = lim -Gb{uj. 0), (3.25) 

a{uj) = -a.(c^,g = 0) = -Gj(u;,0). 
(jj iij 

The DC conductivity can be obtained by further taking the limit 

z % 
(Tdc = lim —Gzz{<^iq = 0) = lim —Gjioj.Q). (3.26) 

Note that all these formulae involve the g = limit. This is a useful feature which 
we will exploit in solving for the hydrodynamic modes from gravity. We will also be 
interested in the thermal conductivity of the charged Dl-brane fluid. The thermal 
conductivity can be evaluated using its relation to the electrical conductivity [17], 
which is given by 

4. Hydrodynamic modes in gravity 

In this section, we study linearised fluctuations of the gravity solution in (2.6) and 
isolate the gauge invariant combinations of fluctuations which correspond to the 
sound mode and the diffusion mode. These we have obtained in the previous section 
using general hydrodynamic considerations. We consider linearised wave like pertur- 
bations of the single charged Dl-brane solution of the form g^i^ — )■ g^i^ + Sg^j,^, -t- 
Afj_ + 6Af^,(j) (f) + S(f) and \E' \1/ + Due to translational invariance along 
the Dl-brane directions, we can assume that all the perturbations can be expanded 
using its Fourier mode as 

%,(t, z, r) = e-'('^*-''^) V(r), M^(t, z, r) = e-'^'^'-'^'^a^ir), (4.1) 
6(f){t, z, r) = e-^('"*-«"V(^), ^Ht, z, r) = e-^('"*-^")^(r). 

We further parameterize the radial dependence of the metric and the gauge pertur- 
bations as 

hu = -4Hu, h,z = cj,H,z, hzz = c\Hzz. a,= ^B„ (4.2) 
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where cx and ct are defined in (2.6). We fix the gauge by imposing dg^.^ = 0, 6 A 
0. The hnearized equations of motion for the perturbations are given by 

3r^H^KH"^ + 3rH [{2H + 1) H + {3H + 1) K] H'^^ + QKr^' 
-AHK {3H + 1) np' + Qr^H'^{H K)B[ + 12 {AH - 1) ^ 

-6{H - 1){H - K)Htt - 8 [QH^ {H + 1) + {H - 1){H - K)} ^ = 0, 
tHH':^ + {bH + 2)H[, + 2r\H -1){H - K)B', = 0, 
3r'^H^KHl[ + 3rH [3H {2H + 1) - {H + 1) K] + QKr^' 

2/2 K 



-AHK {3H + 1) r^' - Qr^H\H - 1){H - K)B[ + 12 2i/ + 1 



+6{H - 1){H - K)Htt - 8 [QH^ {H + 1) - {H - 1){H -K)}ip = 0, 



3 u" 



H'B 



H' 



(4 - H) B', 



H f A 



H^--3^ 



H'.,., — H',. 



8 (H-1) H^ , 

HKB'^ +-{2H{2H + 1)- {5H - 2) K} B'^ + '^H[^ 
+ ^—u{qBt + ujB,) = 0, 



1 



K 



{2H+1) 



np 



6 



2 



r{Hu + H, 



8K {2H + 1) + -^{uj'H - q'K) 



I' 



f = 0, 



(4.3) 
(4.4) 



(4.5) 



(4.6) 



(4.7) 



(4.8) 



r^HK^ + [2H {2H + 1) + K {5H - A)] - -HK{H,, + Hu)' 

16 P 

+At^H\H -1){H- K)B't + — ^ {2H^ + K - H) ^ - A{H - 1){H - K)Ht, 



+ 



[u^H - q^K) 



H 



+ A !^H{2H + 3) -3- {A- H){H - 1) 



K 
H 



e = 0. (4.9) 



Here H and K are defined in (2.7). Equations of motion obtained from the variations 
6g^r and 6Ar lead to the following 4 constraints. 



rH{qKH',, - coHHi) + q{2H + 1) {H - K)Hu - \qK{3H + l)^ + 2g|e 



-2r^H{H -1){H- K){qBt + cuB, 



0, 



TH\qHt. + uH, 



2uji - \uH{3H + l)^ 



H 
K 



{H - K) {2H + 1) {loH,, + 2qHt:) = 0, 



(4.10) 



(4.11) 



3rH'^K {3H + 1) H'^^ + 3rH^ {K + 2H+1) H'^^ + ArH'^K {3H + 1) - QrHK^' 
+Qr^H'\H -1){H - K)B[ + 12 {{2H + 1) H + 2{H - 1){H - K)} ^ 
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-8H {6H\H + 1) + {H -1){H -K)}^- 6H{H - 1){H - K)Hu 

+3^^ (^-q'§Hu + 2uqH,, + uj^H,)j = 0, (4.12) 

r^H^ (^B[ + q^B')j + 2qHt, - UJ - H,, + ^^- = 0. (4.13) 

It can be shown that the constraints (4.10), (4.11), ( 4.12) and ( 4.13) are consistent 
with the dynamical equations of motion (4.3), ( 4.4), ( 4.5), (4.6), ( 4.7), ( 4.8) and 
(4.9). That is, on evolving the constraints using the equations of motion, one does 
not generate new constraints. We have verified that on differentiating the constraints 
with respect to r, one just obtains a linear combination of the dynamical equations 
of motion as well as the constraints. 

Though we have fixed the gauge Sg^r = 0, SAr = 0, there are still residual 
gauge degrees of freedom arising from diffeomorphisms x'^ — )■ x'^ + with e'^ = 
e^{r, oj, g)e~"^*+*5z and U (1) gauge transformations ^ + d^x with x = x{^^ l) 
^-iuit+iqz_ Under diffeomorphism, the metric, the gauge field and the scalars transform 
as 

Qt^u gf,u - V^e,, - Vu^f,, (4.14) 
Af^ — )■ A^ — d^e^Ap — e"dfjAp_, 

where e'^{r,u,q) is determined by the gauge condition Sg^r = 0. The residual U{1) 
gauge transformations on a given Fourier mode of the gauge field act as follows 

At^At- tiox, A,^A, + iqx. (4.15) 

Instead of fixing the gauge completely, it is more convenient to work in variables 
which are invariant under these residual gauge transformations. To do this, we first 
work out the change of the fluctuations under diffeomorphisms explicitly. This is 
given by 

2 

Htt -> Htt - ^(^WQ + n^e^), 

1 2 

Htz -> Ht;, + ^(«we^ - iqet), H^^ H^^ - -^{iqe^ - T^^e,.). 
'^x ^x 

2 L^A L^A 

Bt^ Bt- -^TTH^^ + B,^ B,- iq^^eK 



H' 



(f^if-^er, Y^r, (4.16) 

'^R ^R 

where F's refer to the Christoffel symbols of the single charged Dl-brane solution. 
Similarly under the U{\) transformations, the gauge field fluctuations change as 

Bt^Bt- iiox, B,^B, + iqx. (4.17) 



- 16 - 



From the gauge transformations in (4.16) and ( 4.17), we can show that the following 
are gauge invariant variables both under diffeomorphisms as well as U{1) gauge 
transformations. 

Zp = -q —Hit + 2uqHtz + u) H^^ - — yj, 
H on 

2q 

Gp = qBt + ujB, + 

Sp = 2{1- H)ip + 3^. (4.18) 



where 



V = q\K + 2H + 1) -uj'^{3H + 1) (4.19) 



and k is defined in (2.18). Note that the gauge invariant variables given in (4.18) 
are not unique, in fact any linear combinations of the above variables are also gauge 
invariant. 

After tedious but straightforward manipulations, it can be shown that the dy- 
namical equations and the constraint equations can be used to write down 3 second 
order coupled linear differential equations for the gauge invariant variables Zp,Gp 
and Sp. Before we present these equations, we redefine quantities so that we are 
dealing only with dimensionless variables as follows: 

-| = ti, q = rnq, UJ = r^cu, 

Zp = i^nZp, Gp = rl[Gp, Sp = Sp 

= r%{l + k) with k = ^, L = —. (4.20) 

We also define the expression 

at = q^^-uj\ (4.21) 

In the equations below, for brevity of notation, we continue to refer to the hatted 
dimensionless quantities in terms of their original symbols. The equations for the 
gauge invariant quantities are given below where the prime denotes derivative with 
respect to the dimensionless quantity u. 



-u?{E - 1) + ^u?{E + 1){H- K)]Gp 



Z'p 
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+ (9if^ - 2i7 + - uj'^f - 6H{H - 1){H - - cj^)^ 

-2{H - K)q^{q^ - uj^){2H^ ~ H + 1)] Sp, (4.22) 



G"; - j?t [2g^ir^(l - /f) + oo'H{K{AH - 1) - if(2/J + 1)}] 



^lA^y^'''' - + 1) - -^(// + 1)]^;. - j,^;, 

^x^'^"? + ^(i/ - - K){q\2H + l)(2g^K 
4 



VH^K'^at 



-u\H + K)) - u^atH{3H + 1)}] + —^—{2H + 1){H - K)Zp 
1 .[{l + 3H-6H^)H\q^-uy 



SuH^KVat ' 

- K){2H + {3H + 2)HK + 2{H - 1)7^2} + 3q^HV{H -1){H - K) 
WH{H - K){q\H^ - 6KH - 2) + 2Hat + 2V{H + l)}]Sp, (4.23) 

5^ - ^Mlimt^z'p -^iH-l)iH- K)G'p 

uVat u^at 

-^iH{2H + 1) + KiH - 2)]S'p = 



(// - l)(/f - K)Zp - ^^tI— (i/ - 1)^(// - K){q' - uj')Gp 



1 



u'^H^KVat 



AK 

-{q^ - u^fH^m'' + -7H + 2) 

+q^H{H - K){at{H + l)(4i/2 -9H + 6) + q^{8H^ + -7H + 2)} 
+u^{H - K){q^H{2 + IIH - 13H^ - 8H^) - H{H - l)at{2H^ - 3H - 6) 
+Auj^H{H -1){2H +l)}]Sp. (4.24) 

At present, it seems that there are 3 gauge invariant modes in contrast with the 
2 modes in 1 + 1 hydrodynamics as shown in the previous section. We will show 
subsequently that one of these modes can be decoupled from the rest and consistently 
set to zero and plays no role in determining the transport coefficients. 

4.1 Properties of the fundamental equations 

Though the equations given in (4.22), ( 4.23) and ( 4.24) seem a set of complicated 
coupled differential equations, we will show that for the transport properties of in- 
terest, namely the conductivity and the bulk viscosity can be obtained from them 
using analytical methods. For this purpose, we need to discuss various properties 
relevant to these equations. 
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(i) / = limit 



In this limit, the charged Dl-brane reduces to the uncharged Dl-brane. An 
important check for the system of equations in (4.22), ( 4.23) and ( 4.24) is that they 
decouple and one of the mode reduces to the sound mode studied in [7]. Setting 
/ = 0, we see the parameters which enter these equations reduce to 



ro th, H ^ 1, 
^^g'(3 + /-4A), 



K^f = l 



with A 



(4.25) 



(4.26) 



Note that with these parameters, the mode Zp reduces to the sound mode studied 
in [7], also the definition of Hu is negative of Hu in [7]. Substituting these values of 
the parameters into the fundamental equations for the gauge invariant fluctuations, 
we see that the variable Sp can be consistently set to zero and the equation for Zp 
decouples from Gp. The equation for Zp reduces to 



6 + / 12(1-/) 



rf r(3 + / - 4A) 



(4.27) 

It can be seen that this is the equation for the sound mode obtained in [7] . To write 
the equation for the gauge fluctuation in the / 
it as 



36(1 - ff 
r2/(/ + 3-4A) 



limit, it is convenient to redefine 



Hf 



q 



Bt 



L2 

q6At + uj6A^ + 



Ir'n 



L2 
2qlrl 



2qlrl 



(4.28) 



Thus in the / — )■ limit, the dilaton fluctuation decouples from the gauge invariant 
combination Hp. Substituting Gp in terms of Hp, it can be seen that the sound 
mode Zp decouples from the gauge mode and reduces to 



H'^ 



(A 



4/^ 



q'm, 



0. 



(4.29) 



It can be easily verifled that this is the equation which is obtained by examining the 
gauge field equation 

df^i^e-^'/'/^F'"') = 0, (4.30) 

where the metric and the dilaton background values are that of the uncharged Dl- 
brane. The background gauge field in this case vanishes and the field strength F'^'^ 
is just that of the fiuctuations 6Az and 6At. Thus we have seen that in the / = 
limit, we obtain two modes, the mode Zp corresponds to the sound mode and the 



-19- 



mode Hp corresponds to the charge diffusion mode. This is what is expected for the 
uncharged Dl-brane. The dispersion relation for the quasi-normal mode of Zp was 
obtained in [7] and it is given by 

u; = ±^q-t^q' + .... (4.31) 

Note that we are measuring all quantities in units of th = here. Then identifying 
the sound speed and the bulk viscosity from the above dispersion relation, it was 
seen that 

where s represents entropy density for the uncharged Dl-brane. In this paper, we 
will show that the ratio (/s continues to be ^ for the case of the charged Dl-brane 
also. The quasi-normal mode for the gauge field equation (4.29) is given by 



t3 

Hp = A{1 - u^y'"^"^ { 1 + tu — 



2'"^^ (4.33) 



where A is an arbitrary constant. Note that the above solution satisfies the ingoing 
boundary condition at the horizon u = 1. Imposing the Dirichlet condition at the 
boundary m = 0, we obtain the charge dispersion relation 

uj = -i—q^ + ---. (4.34) 

Here again, we are measuring all quantities in units of th = tq. Using the expression 
for the charge diffusion constant in terms of conductivity given in (3.16), we find the 
conductivity for the Dl-brane system in absence of charge density is given by 



ii. g = limit 



Note that the formula for conductivity as well as the Kubo's formula for bulk 
viscosity involves the g — > limit. It is therefore useful to examine the fundamental 
equations in this limit. The following simplifications occurs in this limit 

V ^ -(3i7 + l)a;l (4.36) 

Examining the equation for the gauge field ( 4.23), we see that it decouples from Zp 
and Sp and it reduces to 
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The equation for Zp and Sp are coupled, they reduce to 

U<-2H-^ %J^X (4,38) 
_{H-l){H-K) {9H'-2H + 1)-6{H-1){H-K) ~ 



{K{H-2) + H{2H+l)} ~ 6{H-1) 
uHK ^ u{3H + l] 



3{H-l){H-K)^ fL^co'H 2H , {H^-l)i3H-2) ^ ^ 

H H 0,-^0/„:r^ . ,X r 'JP; 



M2ir(3/7 + l) [47^2 M2i^ u^H^{3H + l 

where 

5p = w^Sp. (4.39) 
It is now possible to decouple the equations for Zp and Sp by redefining Sp as 

In terms of Sp, the equations in (4.38) reduce to 

^„ {K{H-2) + H{2H + l)} ^, 
uHK 

L^u'^H 2H (if2_ i)(3ij_2) . 



4^2 v?K u^m{3H + l) 
K-2H-1 2{H-l) \ {H-l){H-K) L^u^H 

uK uH{3H + 1) J ^ u^HK{3H +1) ^ 4A'2 ' 



(9i/2 - 2i/ + 1) - 6(i7 - l)(i/ - K) (^^ ^ mi-H) ^^ J ^^^^^^ 



3v?H^K{3H +1) V 3i7 + l 

Note that Zp decouples from the equation for Sp. We can now set Sp consistently 
to zero and study only the decoupled equation for Zp. Simplifying the equation for 
Zp, we obtain 

\ K-2H-l , 2{H-l) \ ^, , L6a;2g ^ 

(g - l)(K(3g - 7) + (3g + l)(2/7 + 1)) 
= «2M(3i7 + 1)2 ^4-42) 

Thus we have shown that setting g = 0, we can obtain two decoupled equations 
(4.37) and (4.42) which correspond to the charge diffusion mode and the sound 
mode. Thus to obtain conductivity and the bulk viscosity of the charged Dl-brane 
fluid, it is sufficient to study the equations (4.37) and (4.42). 

iii. Behaviour at the horizon 



- 21 - 



To obtain the behaviour of the functions Gp and Zp at the horizon, we define 
X = ln(l — u). Then both the equations (4.37) and (4.42) reduce to the oscillator 
equation 

The ratio, (^2k+3)^ obtained due to the behaviour of the coefficient proportional to 
cj^ in both the equations. Thus the behaviour near the horizon is given by 

Gp, (1 - m)=^' 2(2^5)^^ for u^l. (4.44) 

Since classically horizons do not radiate, we need to choose the ingoing boundary 
condition 

(1 - u) *2(2fc+3) ^ (4.45) 

to solve these equations. 

iv. Behaviour at the boundary 

Examining the coefficients of the equation for Gp given in (4.37) for m — 0, the 
boundary, the equation reduces to 

G'; + 2kG'p + —u'^Gp = 0. (4.46) 

Thus the solution for Gp at the boundary, m — 0, admits a Taylor series expansion 
of the form 

Gp^ A{l + 0{u'^)) + Bu{l + 0{u^))+ n^O, (4.47) 

where A and B are integration constants. Similarly examining the coefficients of the 
equation for Zp given in (4.42), we see that, at the boundary, the equation reduces 
to 

- -z; + ^zp = 0. (4.48) 

The above equation admits an expansion of the form 

Zp A{1 + ■■■) + Bu%l + ■■■). (4.49) 

The behaviour at the boundary is necessary to obtain the transport coefficients. In 
fact, the transport coefficients are proportional to the ratio B/A, that is the ratio of 
the normalizable mode by the non-normalizable mode. 

5. Transport coefficients from gravity 

We first summarize the method put forward by [21, 22] to evaluate transport coeffi- 
cients from gravity. 
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1. Let Zk{r) be the gauge invariant variables constructed from the fluctuating 
gravity fields. In general, they satisfy coupled second order linear differential 
equations. We choose linear combination Z{r) such that they satisfy decoupled 
second order linear differential equations. These decoupled gauge invariant 
variables correspond to the hydrodynamic modes of the field theory. 

2. A local solution of the second order differential equations near the horizon 
r = tq will in general be a superposition of incoming and outgoing waves. 
Classically the horizon does not radiate, therefore, we choose the incoming 
wave boundary condition at the horizon. 

3. The solution which obeys incoming wave boundary condition at the horizon 
can be written as a linear combination of two local solutions fi{r) and f2{r) at 
the boundary r — )■ oo as 

Z{r) = Af,{r) + BMr), (5.1) 

where A and B are the connection coefficients of the corresponding differential 
equations. Coefficients A and B depend on the parameters u, q which enter 
the differential equation. Near the boundary, the solution (5.1) admits an 
expansion 

Z(r) = A(1 + ---) + Sr-^(1 + •■■), (5.2) 

where the ellipses denote higher powers of r which are suppressed as r — > oo 
and A > 0. 

4. The action of the quadratic fluctuations can also be organized in terms of the 
gauge invariant observables. Evaluating the action on shell, it reduces to a 
boundary term which is of the form 

S'^) = lim / ^i,Fir, ,)Z'(r)Z(r) + contact te.m=. (5.3) 

r—^oo J 

where the contact terms do not involve derivatives of Z{r) and 

F{r,uj, q) ^ r'^'^^ f{u, q), as, r — > oo. (5.4) 

5. We can now use the fact that Z{r) is a linear combination of the fluctuation 
gravity fields and apply the prescription in [21, 22] to compute the retarded 
correlator for the corresponding operator O in the field theory. We obtain 

B r^+i dZ{r) 



{00) 



R 



A Z dr 



(5.5) 
) finite term 



We have not written an equality but used ~ as we have not yet kept track of 
the proportionality constant which depends on F{r,u,q) in the limit r — oo. 
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We have used the expansion in (5.2) and the property (5.4) to write the last 
relation in (5.5). Note that from the last expression in (5.5), we need to extract 
the finite piece to obtain the ratio B /A. 



6. To apply Kubo's formula, we only need the retarded correlator with g = 0. 
Thus it suffices to evaluate the following ratio to obtain the transport coeffi- 
cients of interest, 

1 dZ{r 



r'-^oo Z dr 



(5.6) 

g=0 



In fact, for the DC conductivity and the bulk viscosity, we need to take a a; — )■ 
limit which is given by 



^ . , 1 dZ(r) 

Re lim ^ ^ 



r^oo,uj-^o icoZ dr 



(5.7) 



q=0, 



Since it is only the ratio 

„ , , 1 dZir) 

at r — )■ oo that determines the retarded correlators, one can determine the differ- 
ential equation satisfied by TZ{r) from the second order ordinary linear differential 
equation satisfied by Z{r). We will see that this is a first order, ordinary but non- 
linear differential equation. The boundary conditions for this differential equation 
are determined from the ingoing boundary conditions satisfied by Z{r) at the hori- 
zon. This equation, in fact, governs the radial evolution of the transport coefficients. 
We will show that for the DC conductivity and for the bulk viscosity, this equation is 
exactly solvable enabling us to determine the analytic expressions for these transport 
coefficients. The fact that the evaluation of transport coefficients can be reduced to 
solving a first order but non-linear differential equation has been observed recently 
for the case of A/" = 4 super- Yang Mills by [23] and has been argued to be true in 
general in [24]. 

The rest of this section is organized as follows: We first show that the radial evo- 
lution of the transport coefficients are determined by first order non-linear ordinary 
differential equations. These equations are exactly solvable for the DC conductiv- 
ity and the bulk viscosity. We then evaluate the effective action to determine the 
proportionality constant relating the ratio TZ in (5.8) to the transport coefficients. 

5.1 Radial evolution of the transport coefficients 



Radial evolution of conductivity 



Let us obtain the equation that governs the radial evolution of conductivity. 
Note that the equation for Gp can be written as 



iHKG-y - ^^"iif"^^ gp + —yCp = 0. (5.9) 



One can now think of this mode as a minimally coupled scalar with a mass term 
proportional to {H — 1). Thus except for the term proportional to {H — 1), it falls 
in the class of equations of motion studied in [25] for which the radial evolution of 
the transport coefficients was easy to obtain in the a; — ?■ limit ^. To remove this 
term from (5.9), we perform the following redefinition 

G'p = ^^^G. (5.10) 
Then the equation for G reduces to 

^ + U(2i/ + i) - ^ + Ti^- ^ = °- ^'-^^^ 

Thus, the redefinition in (5.10) removes the mass term and reduces the equation 
to that of a minimally coupled massless scalar. To obtain the R-charge retarded 
correlator, we need to impose ingoing boundary conditions at the horizon, n = 1 on 
Gp. From the redefinition in (5.10), we see that this translates to ingoing boundary 
condition on G. From the discussion in around (4.45), we see that we have to impose 
the condition 

G ~ (1 -n)"*^(5sw, u^l. (5.12) 
As we have discussed earlier conductivity is proportional to the ratio 

iuGp du CO H{2H + 1) iLoGdu ^ ' ' 

where we have used the redefinition given in (5.10) and also changed the variable 
from r to u. Thus we need to evaluate the ratio at the boundary subject to 

the condition (5.12) at the horizon. Let us define this ratio as 

f - ^ (5 14) 

icuG du 

The boundary condition for this ratio at the horizon, m = 1 is then given by 

^^''■-^2(3 + 2A:)(l -«)+•••' ^^-^^^ 



^See equation (38) in [25]. 
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where the eUipses refer to sub- leading terms at u = 1. The differential equation 
satisfied by the ratio fc can be obtained from the differential equation in (5.11). 
This is given by 

/ 8H^ + 1 2H + 1\ ^ L^H . 
fc + . ^.o^ , 1^ - — ^ fc - ^T^^ + ^^/g = 0. (5.16) 



uH{2H+l) uK ) AK^ 

This is a first order non-linear differential equation which governs the radial evolution 
of conductivity. From this equation, it is easy to obtain the DC conductivity and the 
pole at a; — >■ present in the imaginary part of the conductivity. We first decompose 
the above equation into its real and imaginary parts. 

These equations simplify and decouple in the limit w — ^ 0. This decoupling would 
not have been possible in the original equation for Gp given in (5.9) due to the 
presence of the mass term proportional to {H — 1). But removing this term through 
the re-definition in (5.10) enables us to calculate DC conductivity exactly as follows. 
The solution for fc satisfying the boundary condition (5.15) in the a; — )■ limit is 
given by 

We can now use (5.13) to evaluate the ratio which is proportional to the real part of 
the DC conductivity. This is given by 



Re(7^Gp)„_^o,a;_,o 



L 



3 



Re/c + Re 



2 

L3 {2k + 3f 



k 



ujH{2H+V 



(5.1^ 



2 9VlTk' 

Now from the solution for Im/^ given in (5.17) in the a; — )■ limit, we see that 
Im/c = 0{uj). Thus imaginary part of the conductivity in the w — )■ limit is given 
by 

ImPJc,)„-.„,„„ = Im {^JiJ^)^^^ = (5.19) 

Therefore we see that the imaginary part of the conductivity has a pole at a; — ?■ 
limit which is expected because of the translational invariance of the system. Trans- 
lational invariance implies that there are no-impurities, which in turn implies infinite 
conductivity at a; = by Drude's formula. In fact, using the Kramers-Kronig relation 

lma{u) = --V r^^P^du\ (5.20) 



-26- 



we see that the real part of the conductivity contains a delta function if and only 
if the imaginary part has a pole. Since we have found a pole in the imaginary part 
of the conductivity, it follows that the real part has a delta function singularity at 
a; = 0. Therefore the value for the DC conductivity ^ is valid at w — )■ 0"^. 

As a further check on our analytical manipulations, we have solved the differential 
equation for conductivity given in (5.9) numerically subject to the ingoing boundary 
conditions at the horizon and evaluated the ratio TZcp- For very small values of w, 
we find very good agreement with the formula given in (5.18) and (5.19). This is 
shown in figure. 1 of section 5.2. 



Radial evolution of bulk viscosity 



The bulk viscosity is determined from the equation for Zp given in (4.42) which 
can be written as 

1 d [r^dZp\ 2{H-1) dZp tu^L^H 

_ {H- 1){K{3H - 7) + {3H + 1){2H + 1)) ^ 

- 9umK{3H + 1)2 ^^-^^^ 

where y = v?. Again, we see that the equation resembles a minimally coupled scalar 
equation except for the terms proportional to {H — 1). We can remove these terms 
by the following redefinition for Zp. 

Zp = ^-^Z. (5.22) 

Then the equation for Z reduces to the simple form 

K -2H-1 cu^L^H 

To obtain the retarded two point function of the stress tensor we need to impose 
ingoing boundary condition on Zp at m = 1. Using the redefinition in (5.22), we see 
that this translates to the ingoing boundary condition on G at the horizon. Therefore, 
we need to impose 

Z ^ (1 -M)"'^^(5@f, M->1 (5.24) 

From the earlier discussion, we see that the bulk viscosity is proportional to the real 
part of the following ratio evaluated at the boundary. 



Z" + — Z' + -j^Z = 0. (5.23) 



^Recently [26] has made a proposal for the value of the DC conductivity for conformal systems 
with chemical potential in arbitrary dimensions. We thank Sean Hartnoll for bringing this reference 
to our attention. 
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^ , - k \ ^ f I dZ 

Re ——r- ^ +Re 



Re 



iuj ?>v?{l + ku) J \iuj?)V?Z du 

1 dZ 



iijj?>v? du 



Here, we have used the redefinition of Zp given in (5.22). We are dividing by 
so that we can extract out the ratio B/A in the expansion of Zp near the boundary 
given in (4.49). Note that since the bulk viscosity is proportional to the real part of 
the ratio . l o ^f^, it is determined by the behaviour of Z. Therefore let us define 
the ratio 

/. = ^f- (5.26) 

luZ du 

Using the ingoing boundary condition for Z at the horizon, boundary condition for 
fz at the horizon is given by 

^^'"-^ ^ 2(3 + 2A:)(l-n) + " " " ' ^^'^^^ 

where ellipses refer to sub-leading terms at m = 1. The differential equation satisfied 
by ]z can be obtained from the differential equation for Z in (5.23) and is given by 

fz + fz - ^J^uJ + zufl = 0. (5.28) 

Again separating into the real and imaginary parts we obtain 

Refz + ~^ Re/z - 2a;Im/^Re/z = 0, 

uK 

l^fz + ^'^!^~^ lmfz + c.(Re/| - Im/I) - ^ = 0. (5.29) 
uK 4ii ^ 

The solution of these equations in the a; — )■ limit which obeys the boundary condi- 
tions in (5.27) is given by 

Re/z = Imfz = 0. (5.30) 

It is now easy to obtain the ratio which is proportional to the bulk viscosity. It is 
given by 



1 dZ 



iLuSu"^ du 

fz_ 

'?,u^ 

L^Vl + k 



(5.31) 



Re(7^Zp)U^o = Re 
Re 

3m^' 

(5.32) 
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Again as a further check on our manipulations, we evaluate the ratio TZzp directly 
by solving the differential equation (4.42) numerically subject to ingoing boundary 
conditions at the horizon. We find the result for a; — )■ in very good agreement with 
the expression given in (5.31). This is shown in figure 2 of section 5.2. 

Note that the problem of obtaining the DC conductivity and the bulk viscosity 
has been reduced to solving a first order but non-linear differential equation. This 
equation governs the radial evolution of the ratio which is proportional to the re- 
spective transport coefficient. In the a; — > limit, the solution of these transport 
coefficients were easy to obtain exactly. 

5.2 Comparison with numerical analysis 

In this section, we solve the equations of motion for the charge diffusion and sound 
mode numerically and find the transport coefficients. We will actually find the ratio 
TZgp and Re TZzp which is proportional to the conductivity and the viscosity. Further 
more, we work in a normalization in which = 2 for convenience. Since we have 
analytic expressions for DC value of conductivity as well as viscosity at very small 
u, we can check our numerics with these results. We also know the exact expression 
for conductivity and viscosity in the limit / = and this gives us another check on 
our numerical results. 
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Figure 1: Plots of real (on left) and log of imaginary part of conductivity vs k for the 
single charged case. On left, the different colors red, black, blue, green, magenta and 
purple correspond to w = 10"^", 0.2, 0.4, 0.6, 0.8 and 1.0 respectively. On right, the 
different colors red, blue, orange, gray and pink correspond to cj = lO"^'^, 10~^, 10~^, 10"'^ 
and 10~^ respectively. The dots are the numerical values and the solid lines are curves Im 
a = a is in units of {16ttG^)~^ and oo is in units of 2r'jj/L^. 

In figure (1), we plot real and imaginary parts of conductivity vs k for the single 
charged case. For k = 0, the real part of conductivity approaches 1. This is in 
accord with our analytic calculation for / = case. We compare the k dependence 
obtained numerically for the real part of conductivity for very small u [u = 10~^°) 
with the DC conductivity. We find good agreement between them as the absolute 
value of the difference between numerical and analytically obtained values is at most 
10~^. We expect that the errors in our numerics remain in the same order for all 
other numerical curves, which tell dependence of AC conductivity on k for different 
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values of oj. We don't have analytic expressions for non-trivial u to compare them 
with. We find that conductivity increases monotonically with thought the slope 
decreases as we increase u. Similarly we see that for small our numerical results 
for imaginary part of conductivity fit well with the analytic expression. The absolute 
difference in this case is at most 10^^. For small cj, imaginary part of conductivity 
grows linearly with k. 




Figure 2: Plots of viscosity vs k (on left) and uj for single charged case. On left, the 
different colors red, black, blue, green, magenta and purple correspond to w = lO"^'^, 0.2, 
0.4, 0.6, 0.8 and 1.0 respectively. On right, the different colors red, black, blue, green, 
magenta, purple and orange correspond to k = 0, 0.5, 1.0, 1.5, 2.0, 2.5 and 2.9 respectively. 
C, is in units of r'^/ {IQttG^L'^) and uj is in units of 2r|^/L^. 

In figure (2), we plot viscosity against k and uj for the single charged case. We 
note from the plot that the smallest value of C is at w = = 0. We also compare the 
analytic expression for the viscosity as a function of k for w — i- with the numeric 
plot of C vs /c for UJ = 10~^° (red curve in left plot in figure (2) ). We find the absolute 
difference between analytical and numerical values to be less than 10"'^. From the 
curves, we see that the curve for ( vs k for a given u shifts as a whole as one changes 
oj. From the right plot, we see that the amount of shift increases non- linearly with 

UJ. 
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Figure 3: Plots of real (on left) and log of imaginary part of conductivity vs k for the 
equal charged case. On left, the different colors red, black, blue, green, magenta and 
purple correspond to a; = 10"^", 0.2, 0.4, 0.6, 0.8 and 1.0 respectively. On right, the 
different colors red, blue, orange, gray and pink correspond to oj = 10^^*^, 10~^, 10~*, 10~^ 
and 10~^ respectively. The dots are the numerical values and the solid lines are curves Im 
a = 1^. a is in units of (IGvrGs)"^ and oj is in units of 2r'jj/L^. 
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In figure (3), we plot the real and imaginary part of conductivity against k for the 
equal charged case. Here too, we have analytic expressions for the DC conductivity 
which we compare with the a vs k plot for u = 10~^*^, red curve in the left plot. 
We find a good agreement with the absolute difference between the numeric and 
analytic values being less than 10~^. This bound on error is also same for the plots 
of imaginary conductivity vs k on the right. We observe here that there is little 
change in the curves of cr vs /c as one changes u. The a vs k behaviour here is very 
different from the same in single charged case. Latter, the curves were monotonically 
increasing, but here, conductivity decreases with increasing k. 




k 



Figure 4: Plots of viscosity vs k (on left) and uj for equal charged case. On left, the 
different colors red, black, blue, green, magenta and purple correspond to uj = 10"^'^, 0.2, 
0.4, 0.6, 0.8 and 1.0 respectively. On right, the same correspond to A; = 0, 0.5, 1.0, 1.5, 
2.0 and 2.5 respectively. Orange curve on right correspond to /c = 2.9. is in units of 
rjj / {IdirGsL'^) and uj is in units of Irf^/L^. 

Now we plot the behaviour of viscosity vs k and u in figure (4). Again as in 
the single charged case, the minimum value of viscosity is at /c = w = 0, which is 
same as before and saturates the conjectured lower bound on bulk viscosity. The 
red curve on the left, which stands for viscosity ys k at u = 10"^° is compared to 
analytic value of viscosity obtained in the w — )■ limit. We get the absolute difference 
between analytic and numerical values to be less than 10~^ here, suggesting excellent 
agreement. We find little dependence of viscosity on u, particularly at smaller values 
of A;. 

5.3 Evaluation of the transport coefficients 

In this section, we use the standard prescription of the gauge/gravity correspondence 
to evaluate the retarded two point functions which determines the conductivity and 
bulk viscosity. This will determine the proportionality constant between the ratios 
T^Gp, T^Zp and the transport coefficients. For this, we first need to expand the 
bulk action given in (2.8) along with the Gibbons-Hawking boundary term to second 
order in the fluctuations H^y, B^, ip and ^. In this section, we will not be using the 
dimensionless variables given in (4.20). All quantities in this section will have their 
respective dimensions, whenever needed, we will restore the dimensions of the ratios 
TZgp and TZzp- Using equations of motion and the constraints (Equations (4.3) to 
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(4.13)), we can write the bulk action (2.8) expanded to second order in fluctuations 
as a total derivative in r. 



V K 3 

-^i^B = —r C^HttHtt — HzzHtt — HttHzz + 2HzzHzz ) — -HHtzH^ 

r' 4 2 

-^{H -1){H- K) [qBt + wBz) [r^HKB'^ + 2Htz) - ^K^^' 

ee' + ^M^' + + 4^(2 - i/)e^ + - i)^e 



2i72^^ m^^^ ' ' ' rH^^ 3rH^ 
^ :HuHzz{H{2H + 1) + {4H + 1)K} 



2rH' 

+ ^{Htt + ^..) {K{?,H + l)i7,, + i7(2i7 + K + l)Hu} 

1 4i^' 
+—Htz^{H{2H + 1) - (5iJ + 2)A'} - —{Hu - HzzW 
rK 3r 

+ ^{(7i^ + l)i^tt - (i^ - l)Hzz}^ - ^(Hu + i/..)e (5.33) 

Note that here the prime denotes derivative with respect to r. The Gibbons-Hawking 
term expanded to second order in fluctuations is given by 

SttCs J 
SL^ 4 



-^V-hKe.t = -T,{K{8H + 3) - H{2H + l)}Hl + SrHHtzK 
r K 

-^{Htt - Hzzf{K{l + AH) + H{2H + 1)} 

-2tK{Hu-Hzz){K-H'J. (5.34) 

(2) (2) 

We now combine the S\^^y^ and Sq{^. Using the constraints, we can rewrite it in terms 
of the gauge invariant quantities Zp,Gp and Sp as follows 

Rr.\ / 7r. Rr.\' 

+ contact terms. (5.35) 



ll 


fZp 


Sp\ 


fZp 






V V ~ 


3H 




3H 



where 'contact terms' represent those terms in the action which do not contain any 
derivatives in r and the counter terms which render the complete boundary action 
finite. Next we define a new variable 

S = Sp+ ^^^'-^^ Zp. (5.36) 
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Note that this is also a gauge invariant variable. It has the following useful property 



S 

S ^ as g ^ 0. (5.37) 

where S is defined in (4.40). Thus on taking g — >■ limit, we can consistently set S 
to zero. We can now rewrite the boundary Lagrangian using S as 



2r^at \r^VH J Kr^VH ^ ' y2 

K K 

:rp^SS' - —iZpS' + SZ'p) + contact terms. (5.38) 

To evaluate the transport coefficients using the Kubo's formula in (3.25) and (3.26), 
it is sufficient to look at the boundary Lagrangian at g — ?■ limit. In this limit, we 
can set consistently iS = 0. So the boundary Lagrangian can now be simplified as 

r t^r^ 1 1 

^ = -Tl \ T^T^GpG'p + --^ZpZ'p + contact terms } . (5.39) 
L' {zr^uj'' 8uj^ J 

At g = 0, the expression for Gp reduces to 

Gp = ojB^ = u-^A,. (5.40) 

"0 

Substituting this in (5.39), the boundary action involving the gauge field can be 
written as 

= Y^ijs I dudq{A<^^^)Hu^nG,\u=o- (5.41) 

Here we have converted the derivative in r to derivative in u and used the definition 
of TZgp- A^z^ refers to the boundary value of the gauge field. This field couples to 
the R-current of the Dl-brane theory by the coupling 

^coupling = ^ j d'x{J'Af^ + J^4°)). (5.42) 

Then using the gauge/gravity prescription, we can obtain the retarded Green's func- 
tion of the R-current by 

^2^(2) 

= -J^- (5.43) 

Using this prescription and the boundary action for the gauge field given in (5.41 ), 
we obtain the following expression for the R-current correlator from gravity 

= -T^r7lT3^^^Gp|«=o. (5.44) 
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Finally we can compute the DC conductivity using the Kubo's formula 



a DC = Re ( lim -Gzz{uj, q = 0) 

= ifi n T3 lim Ile7^GpU=0: 
2tI ( {2k + 3)2 



1 {2k + 3f 



(5.45) 



Here, in the third step, we have used the result (5.18) and reinstated the proper 
dimensions for the ratio TZgp which has the dimensions of length. As a check of 
the final answer note that at A; = 0, it reduces to the value evaluated using the 
quasi-normal mode analysis in (4.35). 

Similarly we can determine viscosity using Kubo's formula. At g = 0, the fluc- 
tuation Zp reduces to 



Zp = cu'H,, + 2 ^^^^^V - (5.46) 

Substituting this in (5.39), the boundary action involving quadratic terms in the 
fluctuation H^z is given by 

4i = J dujdq{Hf})HunzAu=o. (5.47) 

Again we have converted the derivative in r to a derivative in u and used the defi- 
nition of TZzp- H^z refers to the boundary value of the fluctuation. The boundary 
fluctuations of the metric couples with the stress tensor of the field theory by the 
following action [27] 

^coupling = \ j dMHi?T'' + Hifr- + 2HS^T'^). (5.48) 

Then using the standard gauge/gravity prescription, the two point function of the 
stress tensor is given by 

^2^(2) 

Gzzzz = -^ .rr ( .,rr , y (5-49) 

Using this prescription and the quadratic action for the metric fluctuation given in 
(5.47), we see the above two point function is given by 

1 3r^ 

Gzzzz = - -,„ ^ -j^^^'^Zp\u=o- (5.50) 
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We now can compute the bulk viscosity using the Kubo's formula 

C = Re ( lim -G-^z,zz{(^, ? = 0) 

\ w-s>0 (jj 



6r 

^ limRe7^z, 



\u=0 



= (5-51) 

Here again, in the third line, we have used the expression for TZzp given in (5.31). 
In the last line, we have written the expression for ( using the definition of entropy 
density for the single charged Dl-brane given in (2.17). Thus we see that the ratio 
of bulk viscosity to entropy density remains l/An when the charge density is turned 
on. 



6. Properties of the transport coefficients 

We first summarize the results of the transport coefficients of the single charged 
Dl-brane. 

_ 1 (2fc + 3)^ 
^ 167rG3 9v/m' ^ ' 

^-16^^- 

In this section, we restrict ourselves to only the DC conductivity except in subsection 
(6.2). Using these two results, we can find three more transport coefficients. The 
charge diffusion constant is related to conductivity by (3.16) and is given by 

L'^ 3 -2k 

The thermal conductivity is also related to the conductivity by (3.27) and is given 
by 

fe + pVcx r\ (2fc + 3)(l + fc) 

Finally the sound diffusion constant can be obtained by (3.12) and is given by 

Ds = , X = ; . (6.4) 

2(e + p) 12r2,ym ^ ^ 

As we have noted earlier, the ratio of bulk viscosity to entropy density is inde- 
pendent of the chemical potential and is given by 

- = ^- (6.5) 



D,= -^^^j==. (6.2) 
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This property also holds for the equal charged Dl-brane solution as shown in ap- 
pendix B. Using the formula for the bulk viscosity (6.1), the thermal conductivity 
in (6.3), the Hawking temperature in (2.17) and the chemical potential in (2.21), we 
can show the following relationship between these quantities is true 

^ = (2-L)^ (6.6) 

This relationship is more striking when we write the chemical potential /i in terms 
of its dimensions. Note that the normalization of the gauge field we have used in 
(2.8) is such that it is dimensionless. This is convenient for the gravity analysis, but 
it is conventional for the gauge field to have dimensions of inverse length. Since the 
chemical potential is basically the value of the gauge field at the horizon (2.21), it 
must have the dimensions of inverse length. Let us therefore restore its dimensions 
by defining 

A=f. (6.7) 
Then the relationship in (6.6) can be written as 

^ = 4vr^- (6.8) 

This relationship is similar to the Wiedemann- Franz law seen between thermal con- 
ductivity and electrical conductivity. A similar relationship between thermal con- 
ductivity and the shear viscosity for the single charged D3 brane was observed by 
[17]. 

6.1 Transport coefRcients at criticality 

In this section, we discuss the reason for this property as well as behaviour of the 
transport coefficients near the boundary of thermodynamic stability k = 3/2. We 
first note that the charge diffusion constant Dc for the single charged Dl-brane 
given in (6.2) vanishes at the boundary of thermodynamic instability. This indicates 
that this mode becomes unstable at k = 3/2 and for this case the thermodynamic 
instability can be studied by examining this mode more carefully. As we will see in 
appendix B, this feature does not hold for the equal charged Dl-brane. It was also 
not seen in the analysis of [17] for the single charged D3-brane. Thus this feature 
seems to be specific for the single charged Dl-brane and it is worth exploring this 
further. 

To determine the critical behaviour of the transport coefficients at the boundary 
of thermodynamic instability, we follow the analysis done by [17]. We first define the 
dimensionless chemical potential m as 

m = ^ = ^ = ^ (6 9) 

27rTH 2nLTH (3 + 2/c) ' ^ ' ' 
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Note that fi/T is the natural variable that occurs in charge current (3.1). We can 
invert the relation in (6.9) to write k as 

1 - 12m2 - ^(1 - 24m2) 

Thus, we can re-express the transport coefficients as 

1 A - 12m2 - Vl - 24m2\ /I - 4m2 + Vl - 24m2 
~ 167rG3 V 72y2m4 J V 1 + m2 

, ttL^T^ /I + 6m2 + Vl - 24m2\ /I - 4m2 - Vl - 24m2 



k = ^ (6.10) 



AGs V 18 / V 8m2 

TiL^T /I + 6m2 + v/r^24m2\ /i _ 4m2 - ^1 - 24m2 



AGs \ 18m2 y V 8m 

1 r. ;T-^/l + 6m2- Vl-24 



1/2 



Dc = -Vl - 24m2 



m 



247rT V ^2(1 + m2) 

1 /5 + Vl - 24m2 



487rr V 1 + m2 



The boundary of thermodynamic stability lies at = | or rric = -4^. Expanding 



2 ■■■c 

the transport coefficients near this point, we see that the Dc exhibits a square root 
branch cut at the critical point. The other transport coefficients are finite at the 
critical point rric, however their first derivatives including that of Dc diverges as 
(rric — m)~^/2_ Thus the critical index is 1/2 which indicates that the system exhibits 
mean field behaviour. A similar behaviour was observed for the shear viscosity and 
conductivity for the single charged D3 branes in [17]. 

From the above expressions for the transport coefficients in (6.11), note that 
that ( and kt are written as T2/(m) and Tg{m) respectively. This demonstrates 
that the system has a hidden 2 + 1 conformal invariance since the entropy density 
is proportional to T2 in 2 + 1 dimensions. Also note that the charge and sound 
diffusivity can be written in the scaling form ^f{m). The conductivity just depends 
on the dimensionless ratio m and assumes the scaling form /(m). From examining 
the scaling form, it is easy to see that as T — oo, keeping the chemical potential /i 
constant, all the expressions for the transport coefficients reduce to the uncharged 
case as expected. Another point worth mentioning is that on expressing Gs,L in 
terms of the Yang-Mills coupling and the rank N, the transport coefficients ( and 
Kt are proportional to N'^/ \/X. If at all this system holographically describes a 1 + 1 
dimensional system seen in nature, the scaling behaviour of the transport coefficients 
seen in (6.11) is a possible test. 

6.2 Behaviour of conductivity 

In figure (5), we plot the conductivity vs quantity 1/m, which is proportional to 
temperature, if chemical potential is held constant. For the single charged case, we 
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Figure 5: Plots of real part of conductivity vs 1/m for single charged (on left) and equal 
charged case. The different colors in the left plot, red, black, blue, green, magenta and 
purple correspond to a; = lO^^'^, 0.2, 0.4 , 0.6, 0.8 and 1 respectively. The colours in the 
right plot, red, black, blue, green and magenta correspond to uj = 10~^°, 0.4 , 0.6, 0.8 and 
1 respectively, a is in units of (167rG3)~^. 

can't go to lower values of 1/m < 1/ rric. We note that for both the single charged case 
and equal charged case, the conductivity saturates to 1, as temperature increases. 
This is expected from our uncharged brane analysis. As m — )■ 0, the behaviour of 
DC conductivity is 

15 

WnGsaDc — 1 + — + ■ ■ ■ for single charged case, 

1 — 24m^ + • ■ ■ for equal charged case. (6-12) 

The low temperature behaviour in equal charged case is cr^j^ ~ m^^ ~ for m^^ —t- 
0. 




Figure 6: Plots of real part of conductivity vs uj for the single charged (on left) and 
equal charged case. The different colours in the left plot, red, black, blue, green, magenta 
and purple correspond to m = 10~^, 0.15, 0.18, 0.19, 0.2, 2\/6 respectively. The different 
colours in the right plot, red, black, blue, green, magenta, purple and orange correspond 
to m =10-^ 0.1, 0.2, 0.3, 0.6, 1 and 100 respectively, o is in units of (IGvrGs) ^ and uj is 
in units of Ir^jL?' . 

In figure (6), we show the dependence of conductivity on frequency for various 
fixed values of m for the single and equal charged case. Here, the behaviour of the 
curves are in contrast with each other in two cases. While for the single charged 
case, we find the curves fit well with the expression ~ a(m) + ^^'^l^i for some uj 



-38- 



independent functions a(m), b{m) and c(m) of m. On the right, we see that the 
conductivity increases for intermediate values of m as w is increased. 

6.3 The relation to the M2-brane theory 

It has been observed that the thermodynamic properties of the near horizon geom- 
etry of M2-branes is similar to that of the Dl-branes [28, 29]. We now recall the 
thermodynamic properties of uncharged M2-branes and compare them to uncharged 
Dl-branes. These properties were obtained from [18]. The near horizon geometry of 
M2-branes is AdS^ times S"^ , let the radius of S*^ be L' and the Newton's constant 
in 4 dimensions be G4. The thermodynamic properties of non-extremal uncharged 
Dl-branes and non-extremal uncharged M2-branes with non-extremal parameter ro 
is given by: 





Dl-branes 


M2-branes 


s 


1 '■o 
4G3 -L* 


1 ^0 

4G4 L'-i 


T 


3 ^0 
27r L3 


3 4 


e 


1 4 


1 rl 




47rG4 


P= - f 







Table 2. Thermodynamics of uncharged Dl-branes and M2-branes. 

From the equation of state p = |, it seems that the non-conformal Dl-brane theory 
behaves as though it is a conformal theory in 2 -|- 1 dimensions. 

This similarity of thermodynamic properties of uncharged Dl-branes and M2- 
branes was also seen to extend to the transport properties. In [7], it was noted that 
the bulk viscosity to entropy density of non-extremal Dl-branes is given by l/An. 
This fact was explained by the observation in [30]. Consider conformal hydrodynam- 
ics of a charged fluid in 2 -|- 1 dimensions ^. The stress tensor and the current are 

^[30] considered the case of the uncharged fluid and obtained a general relation between conformal 
hydrodynamics between 2a dimensions and non-conformal hydrodynamics in d dimensions. The 
case when a = 3/2, d = 2 corresponds to the relation between Dl-branes and M2-branes. In 
general, the relation found in [30] relates conformal hydrodynamics in fractional dimensions to 
non-conformal hydrodynamics in integer dimensions. 
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given by 

f"'' = eu^u^ + Pii]"^ + u^u^) - 2ria''\ (6.13) 
where a, 6 G {0, 1, 2}, rjab is the Minkowski metric in 2 + 1 dimensions and 

aab = PaPbd(^cUd) -^Pabd -U, Pab = Vab + Ua^b- (6.14) 

Let us now dimensionally reduce these equations with the ansatz = (m^, 0) where 
/i G {0, 1} along with the assumption that there is no dependence along the direction 
2 for any thermodynamic variable. Then the non-trivial components of the stress 
tensor and the current can be written as 

f = (e + p)m^m^ + pg^"" - 27]^^"' - 7]P^'''d ■ u, (6. 15) 

3' = 0, 

where 

a^, = Pi:P:di^pU„) - P^,d ■ n = 0. (6.16) 

To show the above expression vanishes, one can explicitly evaluate the components 
or else use the fact that it is a traceless symmetric tensor in 1 + 1 dimensions and is 
orthogonal to the velocity vector u^. Thus the stress tensor and the charge current 
in 1 + 1 dimensions is given by 

where R is the radius of compactification. On comparing this form of the stress 
tensor to that given in (3.1) we see that we can identify 

e = Ri, p = Re, a = Rd, C = Rv- (6.18) 

The entropy density s in 2 + 1 dimensions is related to the entropy density in 1 + 1 
dimensions by 

s = Rs. (6.19) 

From this, we can conclude that for a fluid dynamics in 1 + 1 dimensions, which is 
related by compactification on a circle of radius R to conformal hydrodynamics in 
2 + 1 dimensions, the relation 

p='- (6.20) 
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will continue to hold true due to (6.18). Furthermore, we have 

- = ^. (6.21) 

s s 

Thus the ratio of bulk viscosity to entropy density in 1 + 1 dimensions is identical to 
the ratio of shear viscosity to entropy density of the conformal 2 + 1 hydrodynamics. 

In [30] , it was shown that the equations of gravity fluctuations for the uncharged 
Dl-brane which determine the hydrodynamical transport coefficients is a dimensional 
reduction of the gravity fluctuations of the uncharged M2-brane background. This 
fact and (6.21) explains the reason why the ratio of bulk viscosity to entropy density 
for the Dl-brane is given by l/An. It also explains the fact that speed of sound 
for the Dl-brane theory is same as that of the M2-brane theory. One expects this 
argument to go through for the charged Dl-branes and this is the reason we observe 
that the speed of sound is l/\/2 and the bulk viscosity to entropy density is l/Air. 
As an evidence for this argument, we will now show that the 3 dimensional truncated 
action given in (A. 19) which supports the equal charged Dl-brane solution can be 
obtained by dimensional reduction of the following 4 dimensional action. 

S = / d'x^, {^'+ir2- L"F^''F,}j , (6.22) 

where (74 and i?4 are the 4 dimensional metric and the Ricci curvature respectively. 
G4 in the four dimensional Newton's constant and V is the radius of AdS^. This is 
the action which admits the solution of the equal charged M2-brane. Note that the 
near horizon geometry of the equal charged M2-brane is just a Reissner-Nordstrom 
black hole in AdS^. We address the equal charged case since the single charged M2- 
brane has not been studied in the literature. We now compactify the action in (6.22) 
using the following ansatz 

ds^ = dsf^+i) + e-l'^dy^, (6.23) 

Ay = 0. 

As usual, all fields do not have any dependence on the compact direction y. Substi- 
tuting this ansatz in the action (6.22), we obtain 

To bring the action in the Einstein form, we perform the following re-definition 

~g^u = e3 V- (6-25) 

We then obtain 

S = ^ / d'x^, (i?3 - l{d<Pr + ^et* - L'^e^^F-^F,.) . (6.26) 
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Now comparing (A. 19) and the above action, we see that they are the same on 
identifying 

$ = expf^0V L'=^, A,^^. (6.27) 



X J ' 2' ' ^ ' L' 

This observation indicates that the supergravity fluctuations which determine the 
transport coefficients of the equal charged Dl-brane theory can be obtained by di- 
mensional reduction of the fluctuations which determine the transport coefficient of 
the equal charged M2-brane theory. As a result, the transport coefficients of the 
M2-brane theory is related to that of the Dl-brane theory. 

Finally, we mention that from (6.18) the conductivity of the M2-brane theory 
is related to that of the M2-brane theory. The conductivity of the equal charged 
M2-brane theory has been evaluated in [15] ^ and is given by 

Note that apart from the dimensions set by 6*4, the dependence of the conductivity 
is identical to that of the equal charged Dl-brane theory given in (B.16). In table 1, 
we have compared the transport properties of the equal charged M2-brane and the 
equal charged Dl-brane. 



7. Conclusions 

In this paper, we have studied the transport properties of the 1 + 1 dimensional 
SU{N) gauge theory with 16 supercharges of the Dl-branes at finite chemical poten- 
tial in the framework of the gauge/gravity duality. We evaluated the bulk viscosity, 
electrical conductivity , thermal conductivity, the charge and sound diffusivity for two 
cases. One in which the chemical potential conjugate to one of the U{1) R-charges 
is turned on and another in which equal charges conjugate to all the 4 Cartans of 
the 5*0(8) R-symmetry are turned on. In both the situations, we find that the ra- 
tio of bulk viscosity to the entropy density is independent of the chemical potential 
and is equal to l/47r. We showed that for the single charged Dl-brane theory, the 
charge dissipative mode becomes unstable and the transport properties exhibit criti- 
cal behaviour at the boundary of thermodynamic instability. We also demonstrated 
that the shear viscosity and thermal conductivity satisfy a relationship similar to 
the Wiedemann- Franz law. We have observed that the transport coefficients of the 
Dl-branes theory is same as that of the M2-brane theory apart from an overall nor- 
malization which determines the dimensions and suggests a plausible reason for this 
behaviour. The summary of the transport coefficients obtained in this paper and 
their comparison with the transport coefficients of the M2-brane theory is given in 

^see equation (83) of [15] and identify as k 
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table 1. A technical result of our analysis is the following: we reduced the problem 
of solving the second order differential equation which determines the transport co- 
efficient to a first order non-linear differential equation. This equation governs the 
radial evolution of the transport coefficient. We were able to solve these equations 
analytically for the transport coefficients of interest in this paper. 

A possible extension of this work is to compute the transport coefficients when all 
the 4 chemical potentials corresponding to the 4 Cartans of the 5*0(8) R-symmetry 
are turned on. This would provide a complete knowledge of the transport coefficients 
of the Dl-brane gauge theory. It will also be interesting to understand the thermal 
stability of the full system with all the R-charges turned on. Another direction is to 
understand the connection of the Dl-brane theory with that of the M2-brane theory 
better. This would involve an analysis similar to [8]. We need to show that the 
hydrodynamic fluctuations in gravity which determine the transport for the charged 
M2-brane and Dl-brane are related by compactification. From the point of the view 
of the theories of the M2-branes and Dl-branes, it is interesting to note that unlike 
the presently unknown theory of the M2-branes, the theory of the Dl-brane is a 
regular gauge theory in 1 + 1 dimensions. We have seen that the Dl-brane gauge 
theory provides physical information regarding the M2-brane theory. It is worthwhile 
to explore and utilize this fact to understand the M2-brane theory further. 

1-1-1 relativistic hydrodynamics occurs in the short time description of the plasma 
formed after highly relativistic collisions [31]. The equation of state of this plasma 
does not obey p = e/2, however it will be interesting to see if the transport properties 
of this plasma show the behaviour seen here. Another area where relativistic 1-1-1 hy- 
drodynamics could be important is in carbon nano-tubes and graphene nano-ribbons. 
These materials can be described as a graphene layer rolled up and a graphene layer 
whose linear dimensions is much larger than that of its width respectively [32, 33]. 
These systems are relativistic since they are obtained by a dimensional reduction 
of 2 -|- 1 dimensional graphene which is described by a massless Dirac equation. It 
will be interesting to compare the transport properties of these materials with that 
of the field theory studied here. The system we study here has a gap set by the 
Yang-Mills coupling. Hydrodynamics of other 1-1-1 dimensional systems with a gap 
have been studied in [34, 35] ^. Even though we have analysed only a bosonic sys- 
tem, we can think of it describing a 1+1 dimensional condensed matter system or 
a quasi 1-1-1 dimensional system made up of strongly interacting bosonic quasiparti- 
cles which are themselves made up of elementary electrons, just like Cooper pairs. It 
would be interesting to evaluate an effective Lagrangian for such quasiparticles from 
the action of the gauge theory dual to our gravity system and then compare it to 
effective action for one dimensional effective condensed matter systems like Luttinger 
liquids. A curious observation is that our plots for conductivity vs temperature and 

^We thank Subir Sachdev for bringing these references to our notice. 
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frequency for equal charged case qualitatively looks similar to a system of carbon 
nanotubes-polyepoxy composites [36]. 
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A. Consistent truncation to 3 dimensions 

We first show that the solution (2.6) in 3 dimensions is a consistent truncation of 
the spinning Dl-brane solution in 10 dimensions given in (2.2). For this, we use 
the results of [19] who gave the most general ansatz for the consistent Kaluza-Klein 
reduction of a 10 dimensional solution on the seven sphere The ansatz is as 
follows: 



F(3) = F^ + F^ + F\ (A.i; 



where 



F^ = g~X",' * 'I^Tj, A (/P/i* 
1 



- --^,Tg%,UF'^^AViJ^AVi^\ 

Vn' = dfi' + gA'^fi^, 
VT,, = dT,, + gA^'n, + gA^'Tk,, 

Fi^^ = dA'^ + gA''' A A''^, (A.2) 

and 

/iy = l Ac = T,,/iV t/ = 2T,fcT,fc/iV - AcTii F = det(T,,). (A.3) 



i°See section 5. of [19]. 

^^Note that the sign of here is negative of that in [19], this is a result of a different convention 



for the volume form 63. 
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* is the Hodge dual in the three dimensions. The /x's are defined as follows, 



ui = sin 6', U2 = cos6 simjji, 
u-s = cos 9 cos ipi sin ip2, 1^4 = cos 9 cos ipi cos ■ip2 , 

H2a~l = ^aSinCpa, H2a = J^a COS (pa- (A.4) 

Here, a = 1, ■ ■ ■ 4 and i, j = 1, ■ ■ ■ 8. Then [19] shows that on substituting the above 
ansatz in to the ten dimensional equations of motion, there is a consistent reduction 
to the equations of motion for the three dimensional fields. The equations of motion 
for the three dimensional fields can be derived from the following three-dimensional 
Lagrangian: 

C = R*1- ^dYAdY- ^fr.^ * Vf^u A f^i' A Vfu 

-\y-"'Tr,%' * A Fj') - ^FV4|2T.,.r,, - {fuf} * 1, (A.5) 

where % = Y-^/^Tij. 

Single charged Dl-brane 

We now show that the spinning Dl-brane solution in 10 dimensions given in (2.2) 
can be written in the form given in (A.l). For this, we choose 

H = l + —, Tij=X^i)Sij, 

etzr = 1. (A.6) 
For convenience, we also write down the following 

A ^ H 



I' 2^ L'' 1 



A = 1 + — cos' 9, Y 



r 



.16 ff2 



C = -2;!^ (3+^ 

This results in the metric 
ds^ = H-^^\-fdf + dz^) + hI'^ Qdr^ + r\M9'' + A sin^ ^^0^ ^ ^^^2 ^^^2^ 

r 3 3 

-2H-'^'^^lsin^9dtd(t), {kl 
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which is same as that given by equation (2.2). Note that the exponent of dilaton 
is negative of that given in main text. This is just due to difference in conventions 
between [18] and [19]. The dilaton and the three form are given by 



HI 



1/2 



'(3) 



'L6 



3 + 



2/2 \ 

— cos^ 6 \ dt ^ dz h dr + 2sm.6 cos 6— rdt A dz A d9 



L6 



,r^/ 



-2sin6'cos6'-^ci;2 Ad9 A 



(A.9) 



This also agrees with the expression for the dilaton modulo the sign and the two 
form gauge potential given in (2.2). One can also check that by reading out the three 
dimensional metric by comparing (A.l) to (A. 8) , we obtain the three dimensional 
truncated solution in (2.6). We can now proceed to obtain the three dimensional 
Lagrangian for the Kaluza-Klein ansatz in (A. 6). We first define the scalars Zi = 
Y-^/^Xi = Y-^/^X2 and Z2 = Y-^/^Xj for j 1,2. Then the Lagrangian in (A.5) 
reduces to 



C 



R 



1 



32Y 



d.Yd^Y 



2 



2 \Zl 



d^Z^d^Z^ 



2'^A' 



3 

~Z. 



^d^Z^d^Z^ 



12 



4ri/4 zl 



^F^,F>^'' + —Y'/^Z2{Zi + Z2 



L2 



(A.IO) 



On identifying 



g|0 



(A.ll) 



the above action reduces to the one given by equation (2. 



Equal charged Dl-brane 



We now wish to obtain the truncated 3 dimensional solution as well as the action 
when one turns on equal charges along the 4 Cartans of the 50(8) R-symmetry. We 
start with the 10 dimensional Dl-brane solution with equal spins along the 4 Cartan's. 
This is given by [18]. 



ds' 



Hr 



^3/4, 



-fdt^ + dz"^' 



Hr 



-3/4 J- 3 3 



/l3 



-/ ^ ufdtd(t)i, 



i=l 



A. 



H~^dt+^^lY, 



1=1 




A dz. 
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h = l + 



1/2 
2 ' 



Ho 



/ = 1 



r'^h[dO'^ + cos^ 6'(iV^i + cos^ 6^ cos^ V'irfV'2 + sin^ 



+ cos^ 6 sin^ i'id(f)l + cos^ 6* cos^ ipi sin^ V'2'^03 

+ COS^ 6^ COS^ t/'i COS^ '?/'2'^04]- 



(A.12) 



We will now compare the 10 dimensional solution with the form of the Kaluza-Klein 
ansatz given in (A.l). For this, we first assume that the three dimensional metric is 
of the form 

r . _ ^ J 

(A. 13) 



dsl = Z 



-fdf + dz' + ^dr^ 
hf 



and 



where a^'^ 



Tij = ^5ij, A'^ = a{r)a'^dt, 



(A. 14) 



-cr 



2a,2a-l 



1 and zero otherwise, /ij are given in (A. 4). With this 



ansatz, the gauge field is given as 



gUZ 




dt A dz A dr 



dz A d(j)a A VadUa- 



g2^2 \ 13 \l Zf 

Comparing with field strength of solution in (A.12), we get 



(A.15) 



1 



- - - — 7 = 2^r^/)2 

T^K^^ f g^^rhZH2' g 



(A.16) 



By comparing the metric in (A.l) and the spinning Dl-branes solution (A.12), we 
get 



3/4 
2 5 



A solution to the equations in (A.16) and (A. 17) is given by 



(A.17) 



9 = L-\ 

r'^h ' 
/ = /• 



$-1 = g'^r^h, 
Z = {g'r'h)\ 



(A.18) 
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Now using these equations, the effective 3-dimensional action as given by equation 
(A. 5) reduces to 



2 1 $2 

R d (^d^(^ F F^" + 24 



(A.19) 



with the 3 dimensional solution 



dsl = 




-fdt 


h = 


1 + ku, 


A = 


rll 


udt, 


f = 




$ = 


h ' 





dz' 



Ah^f 



du 



(A.20) 



We have changed the radial variable to n = where th is the radius of the horizon. 



Here, k 



We have divided by th to turn some quantities like L, ro above to 



be dimensionless. The variable / above is same as / in earlier part of the analysis. 
The radius of the seven sphere is L. Parameter ro is related to k and the radius of 
the horizon as 



;i + kYr'^. 



(A.21) 



Note that there is no extra scalar in this case. The general compactification given 
in [19] contains 36 scalars, one singlet under S0(8) and the rest which transforms 
as 35. In this equal charged case, we turn on only the singlet which is the dilaton. 
In the single charged case, one more scalar ^ is turned on and this explicitly breaks 
the 5*0(8) symmetry. 



B. Transport coefficients for the equal charged Dl-brane 

In this part of the appendix, we evaluate the conductivity and bulk viscosity for the 
equal charged Dl-brane. We will be brief here since we have provided the details for 
the single charged Dl-brane in the main text. 

We first provide a table listing the thermodynamic properties of the equal charged 
Dl-brane 
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Hawking Temperature (T") 


f-o (6-2k\ 




i jiiiii tJijy J_yv!'lloluy \0 I 


1 rfrn 


AGs 




1 




Pressure =(-free energy density(/)) 




SttGs L"^ 2 


Charge Density(p) 






Chemical Potential(/i) 


irj^(l+fc) 
L2 



Table 3. Thermodynamic properties of the equal charged Dl-brane 

In evaluating the above thermodynamic quantities, we have used the relation in 
(A.21). Note that the Hawking temperature of this black hole is given by 



T 



6 -2k 



A-kL^th V 1 + A; 

From this expression, we see that the black hole is stable only for 

< 3. 



(B.l) 



(B.2) 



We can also examine the Hessian to see if the equal charged solution admits the 
thermodynamic instability seen in the case of the single charged solution. Using the 
expression of the Hessian given in (2.23), we obtain 



8GlL\k + 3) 
4(1 + A;)2 ■ 



(B.3) 



Since k > 0, the Hessian for this case is always positive and therefore this solution 
does not exhibit the usual thermodynamic instability. Thus the range of the allowed 
values of /c is < /c < 3. This is the same range found for the case of equal charged 
M2-branes [37]. Using the expressions for the thermodynamic variables given in table 
3, we can evaluate the relationship between the charge diffusion constant and the 
conductivity from the formula in (3.14). It is given by 



3{k + 3) 
'2rl{3-kf 



(B.4) 



Hydro dynamic modes from gravity 



To obtain the two hydrodynamic modes of the charged fluid from gravity, we 
analyze linearized wave like perturbations in the background of the equal charged 
Dl brane solution given in (A. 20). It is a solution of the action given in (A. 19). The 
perturbations are defined as follows: 

g^^ ^ gl{l + if,,), A ^ A? + 

$->$° + i:V (B.5) 

where the superscript '0' refers to the background values. Due to translational in- 
variance along the t and the z directions, we can assume that the dependence of the 
perturbations along these directions is of the form as ~ exp[-p-(— cut + qz)]. Note 
that here we will be using the dimensionless variables defined in (4.20). 

To write the gauge invariant modes, we first introduce the following functions 

V = q\4h^ - r^y) - 4uj^h\ 

a = q^f-uj'^. (B.6) 

The two gauge invariant variables which are invariant both under diffeomorphism as 
well as f/(l) gauge transformations are given by 

Zp = -q^fHtt + 2a;giite + u'^H^^ + y^v^, 

Gp = qBt + ujB, + q<^. (B.7) 

These gauge invariant variables satisfy the following equations of motion 

VafuZp - 2h^{q^ - u'^){q^{h'^ + 5h - 12) - 2huj^}{l - f)Zp 

+q^h\l - f f{q\2h'' + 2/i - 8) - w''{h'' + 6/i - 8) + q^fh{2>h - ^)]Z'p 

+8ih - 2)h\q^ - u^fZ'p - 2qkrl'^^G'p 
= -8A:gr^^ ^'^y ^\ ?,q'f + Auj' - qW{A + 3f-h + fh)]Gp 

+'^m'-^'f-q'{i-f){h+m'-u') 

+q^h{l - ff{3u~'h'f{h - 4) + q')]Zp, (B.8) 



G"; + ^{g'(2 + h-hf)- 2lo^}Z'p + 



2k h){l-f) 

-T—[q - uj -\ 7-7 

ha anju 



G p 



+ J^.Gp. (B.9) 
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Note that these equations decouple in the q ^ hmit. 



Conductivity for the equal charged case 

In the limit q — 0, the equation for the gauge invariant current mode is 

We can remove the coefficient of Gp proportional to (1 — /) by the following redefi- 
nition 

4 - h 

Gp = —G. (B.ll) 

As wc have seen for the case of the single charged solution conductivity is essentially 
determined by the ratio 

^ 1 dGp ,^ , 

luG p au 

which in turn is determined by the ratio g = Ingoing boundary conditions 

at the horizon for Gp corresponds to the following boundary condition on g 

\im q(u) — — —- r. (B.13) 

^ {3 - k){l + k){l - u) ^ ' 

The equation of motion satisfied by g is given by 
The solution for the real part of gf in a; — )■ limit is 

Just as in the previous case for the single charged, the real part of DC conductivity 
here is proportional to the value of ReT^-G^ at the boundary {u — 0), which is given 
by 

Note that here we have reinstated the factors of which we have absorbed in 

defining cj. The imaginary part of KeTZcp is given by 



L3 

Im7^Gp = yl"^ 



d 4 - h 
du icuh 



= (B.17, 



In Figure 3. we have compared these expressions with that determined by numerically 
solving the equation for Gp. We find that they agree to less than one part in 10~^. 
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The DC conductivity is related to ReGp by a proportionality constant which can 
be determined from the boundary effective action as was done in the single charged 
case in section 5.3. This results in 

^ ' (B.18) 



IGTrGs 9(1 + kf 



Bulk viscosity for the equal chcirged case 

In this section, we calculate the bulk viscosity for the equal charged case and 
show that the ratio ^ is constant. In the g = limit, the equation for the sound 
mode decouples from the current mode. It turns out to be 

Z'i>{u) + ^-^^Z'Au) - ^-^Z'M + ^,Zp{u) = 0. (B.19) 



Let us define the ratio 



9=^y (B.20) 
iujZp[u) 



Since the sound mode satisfies ingoing boundary condition at the horizon, the func- 
tion g should satisfy 

limg=7 , , r. (B.21) 

The appropriate solution for g( in a; ^ limit is 

Reg = {l + kf—. (B.22) 

The bulk viscosity is proportional to the real part of the following ratio evaluated at 
the horizon: 

We can evaluate this from the expression for Reg in the a; ^ limit which results 
in 

Here we have reinstated the factor of which we have absorbed in the definition 
of uu. We have verified that the above expression using the numerical solution for 
the equation for Zp to one part in 10~^. This is shown in figure 4. Evaluating the 
proportionality constant relating the bulk viscosity to the ratio RelZzp, we obtain 
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The entropy density of the equal charged solution (A. 20) is given by 



rlru 



(B.26) 



Using (A. 21) and the expression of the bulk viscosity (B.25), we get the ratio 



c 

s 



1 

An' 



(B.27) 



We can now evaluate the thermal conductivity of this solution using (3.27), this 
results in 

rl {l + k){'i-k) 



8LG. 



k 



(B.28) 



It can also be verified that this system also satisfies the Wiedemann-Franz like be- 
haviour. 



(B.29) 



The remaining transport coefficients Dc and Dg which are related to the conductivity 
and the bulk viscosity can be evaluated and are listed in table 1. In the end, we 
mention that we have verified that the transport coefficients of the equal charged 
solution does not exhibit the critical behaviour seen in the case of the single charged 
solution in the domain of < /c < 3. This is consistent with the fact that the Hessian 
does not show any sign of thermodynamic instability. When written in terms of 



m 



Vk 



2tx LT 3-k' 



the various transport coefficients are 



~ag7 



1 + 2^1 + 12m2 
3(1 + 16m2) 

1 + VI + 12m2 
6 



l + 6m2 + Vl + 12m2 
18m2 



Vl + 12m2 (1 + 24m2 - Vl + 12m2) 



247rT 



m2(l + 16m2) 



247rT 



1 + 2V1 + 12m2~ 
1 + 16m2 



(B.30) 
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